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ABSTRACT 

We present constraints on the number of Galactic magnetars, which we have established by searching 
for sources with periodic variability in 506 archival Chandra observations and 441 archival XMM- 
Newton observations of the Galactic plane (|6|<5°). Our search revealed four sources with periodic 
variability on time scales of 200-5000 s, all of which are probably accreting white dwarfs. We identify 
7 of 12 known Galactic magnetars, but find no new examples with periods between 5 and 20 s. 
We convert this non-detection into limits on the total number of Galactic magnetars by computing 
the fraction of the young Galactic stellar population that was included in our survey. We find that 
easily-detectable magnetars, modeled after persistent anomalous X-ray pulsars (e.g., either with X-ray 
luminosities Lx=3x 10 33 erg s _1 [0.5-10.0 keV] and fractional rms amplitudes A rms =70%, or Lx=10 35 
erg s -1 and A rms = 



12%), could have been identified in «5% of the Galactic spiral arms by mass. If we 



assume that there are 3 previously-known examples within our random survey, then there are 59^32 
in the Galaxy. Barely-detectable magnetars (Lx = 3 x 10 33 erg s _1 [0.5-10.0 keV] and ^4 rms =15%) 
could have been identified throughout wO.4% of the spiral arms. The lack of new examples implies 
that <540 exist in the Galaxy (90% confidence). Similar constraints are found by considering the 
detectability of transient magnetars in outburst by current and past X-ray missions. For assumed 
lifetimes of 10 4 yr, we find that the birth rate of magnetars could range between 0.003 and 0.06 yr -1 . 
Therefore, the birth rate of magnetars is at least 10% of that for normal radio pulsars. The magnetar 
birth rate could exceed that of radio pulsars, unless the lifetimes of transient magnetars are >10 5 yr. 
Obtaining better constraints will require wide-field X-ray or radio searches for transient X-ray pulsars 
similar to XTE J1810-197, AX J1845.0-0250, CXOU J164710. 2-455216, and IE 1547.0-5408. 
Subject headings: stars, neutron — stars, statistics — X-rays, stars 



1. INTRODUCTION 

Astronomers have only started to appreciate the diver- 
sity of the properties of neutron stars that are produced 
when a massive star collapses and explodes (jPopov et al.l 
2006). The list of different manifestations of neutron 
stars now includes: radio pulsars that are powered 
by t he rotation of their 10 8 — 10 13 G magnetic fields 
(e.g..(Lorimer et aL[l2006D; accreting X-ray pu lsars (e.g., 
iBildsten et alJll997Uwii nands fc van der Klislll998l) and 
therm onuclear bursters (e.g., IStrohmaver fc Bildstenl 
12006( 1 that are accreting matter from binary com- 
panions; magnetars that are powered by the de- 
cay of their >10 14 G fields (which historically have 
been categorized as either ano malous X-ray pulsars 
or so ft gamma repeaters; e.g., IWoods fc Thompson! 
2006); intermitt ently-detectable "Rotatin g RAdio Tran- 
sients" (RRATs: iMcLaughlin et~aTll2006[) ; isolated, cool- 
ing neutron stars t hat shine pr imarily in soft X-rays 
(jWalter et alj|1996t lHaberlll2007t ). and central compact 
objects that are seen as poi nt sources near the cen- 
ters of supernova remnants jC hakrabartv et al] 120011 : 
iSeward et all [20031 : IPavlov. Sanwal. fc Teteri 120041 ). Un- 
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derstanding the properties of these compact objects and 
their birth rates provides important constraints on the 
late-time evolution of massive stars, and on the pro- 
cesses that occur during stellar collapse. The relation- 
ships among the different classes of compact object could 
reveal how their magnetic fields decay and their interiors 
cool. 

In this paper, we present a search for magnetars. This 
search is timely for two reasons. First, recent evidence 
suggests that magnetars are the products of unusually 
massive progenitors. Three magnetars have been found 
to be in clusters of massive, young stars dFuchs et al.1 
I1999L IVrba et alJl2000t lEikenberrv et alJ 12004( 1. and the 
turn-off masses of two of these clusters imply that the 
progenit ors to the neutron stars were very m assive, 30- 
40 M Q (jFiger et alJl2005t iMuno et al.ll2006h . A fourth 
magnetar has been associated with a bubble of neutral 
hydrogen that was probably blown b y the wind of a >30 
Mq progenitor ((Gaensler et al.1 12005) . This suggests that 
massive stars may be more likely to produce magne- 
tars, whereas ordinary radio pulsars are generally pre- 
sume d to be left by lo wer-mass, 8-20 M© progenitors 
(e.g.. iHeger et all 12003). Giv en that less m assive stars 
are much more common (e.g.. lKroupall2002Tl . if massive 
stars produce magnetars, one would expect that their 
birth ra tes should be much lower than those of radio 
pulsars ((Gaensler et alJl2005h . 

Second, there is significant debate about how the 
strong magnetic fields that characterize magnetars are 
produced. The original hypothesis is that magne- 
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tars are born with millisecond periods, and that the 
strong fields are produced by a dynamo in the rapidly- 
rotating proto-neu tron star ((Thompson fe Duncaniri993t 
IHeger et al.ll2005[ ). However, observations of supernova 
remnants associated with magnetars rule out the ex- 
pected input of energy from neutron stars with ini- 
tial spin periods <3 ms (|Vink fe Kuiperl [2006D . At the 
same time, the discovery of a few OB st ars with 10 3 
G surface fields (|Donati et al.l [20021 [2006aH bl) has moti- 
vated the alternative hypothesis that magnetar-strength 
fields are primordial, having been amplified only by th e 
collapse of the core (Fcrrario & Wickramasinghc 2006). 
This second hypothesis is attractive because it makes a 
straightforward prediction, that the birth rate of mag- 
netars should be equal to that of highly-magnetized 
OB stars. However, it cannot yet explain why some 
highly-magnetized neutron stars are radio pulsars instea d 
of magnetars (e.g.. iPivovaroff. Kaspi. fe Gotthelfl 12000) . 
This discrepancy is one of the main points in favor of 
the a-Q dynamo process that would act in a rapidly- 
rotating proto-neutron star, because it could produce 
^,10 15 G internal fields that would powe r the magnetars 
ijThomspon. Lvutikov. fc Kulkarnill2002t ). 

Here, we report the results of our search for new Galac- 
tic magnetars in archival observations taken with the 
Chandra X-ray Observatory and the XMM-Newton Ob- 
servatory. Previous searches for periodic sources with 
these observatories have been lim ited to small fields , 
such as the Small Magellanic Cloud (|Macomb et al.ll2003l 
lEdge et "ail l2004f) and the central 20 pc of the Galaxy 
(|Muno et alj [2003) . Our search incorporates observa- 
tions through the entire Galactic plane. Although we 
find four new sources with significant periodic signals, 
none are likely to be magnetars. Therefore, we use the 
survey to place limits on the number of active magnetars 
in the Galaxy, and discuss the implications for their birth 
rates. 

2. OBSERVATIONS 

We attempted to identify new magnetars by search- 
ing for sources with periodic X-ray pulsations in archival 
Chandra and XMM-Newton observations of the Galactic 
plane. We included observations with a Galactic latitude 
|6|<5° that were public as of January 2007. For both ob- 
servatories, we required that some of the data were taken 
in an imaging mode. For Chandra, we further rejected 
observations taken with the gratings in place. For XMM- 
Newton, we rejected all observations shorter than 10 ks, 
because they were often affected by background flares 
throughout their entire duration. In total, we searched 
506 Chandra observations, and 441 XMM-Newton obser- 
vations. The coverage on the sky is illustrated in Fig- 
ure [TJ Within a Galactic longitude of |Z|<10° and a lat- 
itude |6|<5°, which encloses about half of the mass of 
the Galaxy (e.g., using the model in Launhardt, Zylka, 
& Mezger 2002), these observations cover about 6% of 
the sky. Closer to the Galactic plane, within |Z|<10° and 
|&|<0.5°, the observations cover 25% of the sky. 

The raw event lists and calibration data were down- 
loaded from the High Energy Astrophysics Science 
Archive Research Center. 6 We processed the event lists 
produced by each observatory in the standard manner, in 

6 http://heasarc.gsfc.nasa.gov/ 



order to extract events for sources that we could search 
for periodic variability. 

2.1. Chandra Data Preparation 

The Chandra data were processed using CIAO version 
3.4. As we were only interested in the arrival times 
of events, we did not apply the latest calibration. We 
started with the default level 2 event files, and removed 
any time intervals during which flares from the parti- 
cle background caused the total event rate from the de- 
tectors to increase by more than 2a above the mean 
event rate. We then generated exposure maps at a fidu- 
cial energy of 1.5 keV (i.e., the peak of the detector 
effective area), and binned the event lists to produce 
images in the 0.5-8.0 keV band. These were used to 
search for point sou rces using the routine wavdetect 
([Freeman et al.ll2002t ). For computational efficiency, we 
searched a series of three images using sequences of 
wavelet scales that increased by a factor of \J~2: a cen- 
tral, un-binned image of 8.5' by 8.5' searched over scales 
of 1-4 pixels, an image binned by a factor of two to cover 
17' by 17' searched over scales of 1-8 pixels, and an im- 
age binned by a factor of four to cover the entire field 
searched over scales of 1-16 pixels. The source lists from 
each image were combined for each observation, favoring 
the positions derived from the images with the highest 
spatial resolution for sources detected at multiple scales. 
We did not attempt to discriminate real sources from 
detector artifacts, such as the boundaries of the individ- 
ual charge-coupled devices (CCDs), at this stage of the 
algorithm. 

We then extracted events from each source, using a 
radius defined to enclose 90% of the point spread function 
(PSF) at 4.5 keV. In this relatively higher energy 

was used so that we would not exclude photons from 
sources with hard spectra. The radius was dependent on 
the offset from the aim point in arcminutes (6) , and was 
parameterized based on simulations with the CIAO tool 
mkpsf as 

r = 3.27 + 0.34261 + 0.0206» 2 + 0.0196* 3 , (1) 

where r is the radius in pixels (0'.'492). We corrected the 
arrival times of the events to the Solar System barycenter 
using the tool axbary. 

2.2. XMM Data Preparation 

The XMM-Newton event lists from the European Pho- 
ton Imaging Camera (EPIC) were analyzed using version 
7.0 of the Science Analysis Software, 7 CIAO 3.4, and the 
High Energy Astrophysics Software version 6. 3.0. 8 We 
examined the lists from each active imaging detector sep- 
arately, starting with the files from the archive. For most 
observations, the same field was observed with three in- 
dependent cameras: one pn CCD array and two metal- 
oxide-silicon [MOS] CCD arrays. We removed time in- 
tervals during which particle events caused the event rate 
from the detector to flare more than two standard devi- 
ations above the mean rate. This selection was generally 
successful at removing particle flares, but in many cases 
the entire observation for a given camera was affected by 



7 http://xmm.vilspa.esa.es/external/xmm_user_support 
/document at ion/sas_wsg/USG/USG . html 

8 http: //heasarc .gsf c .nasa. gov/docs/sof tware/lheasof t/ 
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Fig. 1. — The locations of archival observations used for this survey. Black squares are Chandra observations, schematically represented 
with boxes that have an area equal to that of the ACIS-I. Grey rectangles are XMM-Newton observations, schematically represented with 
boxes that have an area equal to that of one EPIC MOS. 



flares, which rendered this automatic algorithm useless. 
We discarded data from individual cameras (usually the 
pn) for observations too badly affected by flares. 

We then created images of the 0.2-f2 keV events, 
binned to 4" resolution. The standard data selection was 
applied to make the images, in order to remove events 
near the edges of the detector chips and bad pixels, and to 
reject events that were likely to be cosmic rays (pattern 
>4 for the pn and >12 for the MOS). We then generated 
matching exposure maps. We searched for point sources 
using the routine ewavelet, separately for each observa- 



tion. At this stage of the algorithm, we did not attempt 
to discriminate real sources from detector artifacts, nor 
did we attempt to verify that sources detected in one 
camera were present in the others. We extracted event 
lists from individual sources using the radius produced 
by ewavelet, which generally was «15" and enclosed 
«50% of the PSF. The arrival times of the photons were 
corrected to the solar system barycenter using the tool 
barycen. 

2.3. Identifying Candidate Signals 
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For sources with at least 100 total events from ei- 
ther observatory (including background), we computed 
Fourier periodograms to search for periodic signals. The 
range of periods searched for both instruments was de- 
signed to encompass t hose of known magnetars w ith X- 
ray pulsations, 5-12 s (jWoods fc Thompsonll2006[ ) . How- 
ever, we not e that after our analysis was completed, 
ICamilo et all (|2007h announced the discovery at radio 
wavelengths of 2 s pulsations from IE 1547.0-5408. Pe- 
riods this short would generally only be identifiable in 
XMM-Newton EPIC-pn data with our search. 

Fo r the Chandra data , we used the Rayleigh statistic 
(Z\: iBucceri et al1l!983f ). We searched for signals with 
frequencies between 1.5 times the Nyquist frequency (i.e., 
1.5 x l/2£bin, where ibin is the interval at which the data 
were read out) and 10% of the inverse of the total time 
interval of the observation (i.e., 0.1/T cxp ) to avoid red 
noise. We used a frequency step corresponding to the 
inverse of the total time interval (1/T cxp ). For most ob- 
servations, the data were read out every ibin=3.1 s, so the 
highest frequencies searched correspond to periods of 4.1 
s. Observations typically lasted between 1.2 ks and 120 
ks, so we could identify signals with periods at the upper 
range of at least 120 s, and in some cases 12,000 s. 

For the XMM-Newton data, the Rayleigh Statistic was 
computationally inefficient to compute for pn observa- 
tions of bright sources, so we computed discrete fast 
Fourier transforms. The data were padded so that the 
number of points in the transform was a power of 2, so 
the frequency resolution was generally finer than 1/T exp . 
The maximum frequency considered was the Nyquist fre- 
quency for the data. The pn data was taken with a time 
resolution of at least 73.4 ms, providing sensitivity to 
periods as short as 0.15 s. The MOS data was taken 
with a time resolution of at least 2.4 s, and so our search 
was sensitive to periods as short as 4.8 s. The lowest 
frequency considered was 0.1/T oxp . 

The confidence with which we could exclude that any 
given signal was produced by random noise depended 
upon the number of trial signals examined, Atrial- T ne 
typical Chandra observation lasted 20 ks, and contained 
3 sources with >100 counts, and had a time resolution of 
3.1 s, so that there were iV tr i a i~20,000 trial frequencies 
in the ensemble of periodograms from each observation. 
The typical XMM-Newton observation also lasted 20 ks, 
and contained ~25 sources per camera with >100 counts. 
For each MOS detector, with a time resolution of 2.4 s, 
the typical set of periodograms contained At r i a i«200,000 
trial periods per observation. For the pn detector, with 
a time resolution of 73.4 ms, the typical set of peri- 
odograms contained iVt r i a i«1.8 x 10 s trial periods. In 
total, we searched 1.7 x 10 11 trial periods, the vast ma- 
jority of which were from the high-time-resolution data 
taken with the EPIC pn. 

The powers produced by random noise in a peri- 
odogram in which n measurements have been averaged 
are distributed as a chi- squared function wit h 2n degrees 
of freedom. Following (Ra nsom et alJl2~002h . we refer to 
the measured Fourier powers with n—1 as P mC as, and nor- 
malize them so that the mean power produced by white 
noise is 1. The chance probability that noise would pro- 
duce a signal larger than P mea s can be determined from 



an exponential distribution: 

prob = 1 - (1 - e^ p — a jVt r iaie- p — , (2) 

where the appro ximation is valid for P meas 3>l 
(R ansom et alJ l2002t ). Given the large number of 
trials for our entire search, a signal that had a <0.1% 
chance of resulting from noise must have a power 
-fmeas>32.8. Such a signal would be detected at a 
confidence level equivalent to 3cr over the entire search, 
or 8ct in a single trial. 

However, a signal could also be considered significant 
if it was detected with a lower power in multiple observa- 
tions, as one might hope would occur given that XMM- 
Newton has three separate cameras that observe the same 
patch of the sky. For instance, given two signals P meas ,i 
and P meas ,2 at the same frequency, the chance probabil- 
ity that their sum exceeds some value is a chi-squared 
distribution with 4 degrees of freedom: 

prob « iV tria l(l + P m e as ,l+^meas,2)e~ (P — 1+P — 2) (3) 

(see iRansom et alJ l2002t for the general form for sum- 
ming an arbitrary number of signals). If we take 
Pmeas, 1 = Pmeas, 2, for example, a signal has a <0.1% 
of resulting from noise if it appeared with P m cas,i>17 in 
both observations. In principle, one could devise an algo- 
rithm that searched through all of the periodograms from 
the same source, and sum the powers at each frequency 
to search for signals that repeat in the data. In prac- 
tice, however, the periodograms were not all computed 
with the same frequency resolution, which makes such an 
effort difficult. Moreover, when considering observations 
separated in time, one also has to be concerned that some 
candidate signals drifted in frequency, either because of 
the spin-down of an isolated pu lsar, or Doppler shifts for 
a pulsar in a binary (see, e.g.. IVaughan et aLlll994L for 
further discussion). 

Therefore, we have adopted a simplified approach in 
examining candidate signals, by recording all signals with 
powers with less than a 1% chance of resulting from noise 
in a search of a single source. For Chandra ACIS, the 
threshold power is generally P mca s>17. For the XMM- 
Newton MOS, the threshold power is typically P moas >19, 
whereas for the pn the power is P mcas >26. Any candidate 
signal was inspected to determine its significance. 

With these search criteria, our results up to this point 
were dominated by signals that are non-periodic noise or 
detector artifacts. In both Chandra and XMM-Newton 
data, we detected low-frequency noise from astrophysi- 
cal flares in the count rates of individual sources such 
as pre-main-sequence stars, and from background flares 
that our filtering algorithm failed to remove. From Chan- 
dra ACIS, we detected signals from sources that fell near 
chip boundaries, at the harmonics and beat periods of 
the satellite dither (700s in the x-direction, 1000s in the 
y). For the XMM-Newton EPIC, particularly the pn, we 
found signals with a range of periods that appeared to be 
related to hot columns and chip boundaries, particularly 
in observations with high particle background. We are 
not certain of the origin of these signals from the EPIC. 
The spurious signals introduced by the detector generally 
shared the feature that they could be found in multiple 
sources at the exact same frequency during an obser- 
vation. Therefore, we have removed from consideration 
any signals that appeared in two or more sources on the 
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same detector in the same observation. After removing 
such signals, we found 358 sources with candidate signals 
in the Chandra observations, and 1380 sources (some of 
which are duplicates) with signals from XMM-Newton. 

These signals still turned out to be dominated by 
low-frequency noise and detector artifacts, which could 
be quickly determined by visually inspecting the power 
spectrum. Therefore, for the final step, we scrutinized 
«1700 power spectra by eye to remove the remaining ex- 
amples that were clearly noisy, and to remove sources 
that appeared to be detector artifacts. We defined a 
signal as significant if it had a power P meas >32.8 in a 
single observation or had a power larger than the single- 
observation threshold in two or more observations. We 
found a few sources with significant periodic signals that 
we could not attribute to noise or detector artifacts. 
We describe the previously-known and new sources sep- 
arately below. 

2.3.1. Known Sources 

Most of the signals were from previously-identified 
pulsars. These confirmed that our algorithm 
worked as intended. In our Chandra observations, 
we identified the high-mass X-ray binary (HMXB ) 
4U 1145-619 (|White et alJ ll978t iRutledee et al.ll2005h . 
CXOU J164710.2-455216 (|Muno et all 120061) . and two 
sources toward the Galactic center (CXOC J174532.3- 
29025 1 and GXOGC J174532. 7-290552; iMuno et all 
2003). In the XMM-Newton observations, we iden- 
tified the pulsars Gemin ga (jHalpern fc Holtl 119921 : 
I.Tackson fc Hafoernl 12001 and PSR J1513-5908 
([Seward fc Harndenl 119821), the HMXB Set X-l 
(|Kovama et al.l Il99it Kaplan et all I2007D. and the 
magn e tars IE 1048 . 1-5937 ([Seward, Charles, fc Smald 
[19861 iTieneo et alJ I2005T), 1RX S J170849.0-400910 
dSugizaki et alJ 119971 : iRea et all I2005T). SGR 1806 - 
20 ([Murakami et al.l Il994t jMereghetti et all [20051) , 
SGR 1900+14 ([Hurlev et aD Il999t iMereehetti et at 
' 2006al), XTE J1810-197 inoutburst ([Ibrahim eTal 
2004V jHalplrn fc Gotthell 120051 ). and IE 2259+586 



([Fahlman fc Gregorvlll981[ ). The list above includes 7 of 
the 12 confirmed magnetars in the Galaxy. 

It is notable, however, that several magnetars were 
not detected in our search, despite being the targets of 
archival Chandra and XMM-Newton observations. Here, 
we summarize the difficulties encountered identifying 
several examples: 

IE 2259+586 was not identified with Chandra because it 
saturated the detector during an imaging observation. 
IE 1048.1-5937 and 1RXS J170849. 0-400910 were not 
identified with Chandra because they were only observed 
with the gratings in place. These cases are not a serious 
concern, because such bright sources are rare, and so 
almost never are found serendipitously in the fields of 
Chandra and XMM-Newton. 

SGR 1806-20 and SGR 1900+14 were not identified by 
Chandra while in quiescence. Although their signals were 
present in the data, their powers were below our search 
threshold. All of the above sources were identified with 
XMM-Newton. 

4U 0142+61 was not identified with either Chandra or 
XMM-Newton. This is partly because the source had a 
small pulse fraction (4% rms) , but also because of the de- 



tector modes with which the source was observed. With 
Chandra the source either saturated the detector during 
imaging observations, or was observed with the gratings 
in place. With XMM-Newton, only the MOS2 camera 
was active, and the magnetar only produced a signal 
above the single-observation threshold (P meas >18) in one 
of the two observations. That signal (P moas =24.7) was 
below the threshold for our entire search, and so can not 
be considered a detection as part of our blind search. 
However, had that source been observed with the pn ac- 
tive, we would have identified it. 

IE 1841-045 also was not identified with Chandra or 
XMM-Newton. With a fractional rms amplitude of 13% 
rms, it did not produce a significant signal in the Chandra 
data. The XMM-Newton observations of this source were 
too short (<10 ks) to be included in our search. A longer 
XMM-Newton observation would almost certainly have 
identified this source. 

SGR 1627-41, AX J1845. 0-0258, and XTE J1810-197 
in quiescence were all too faint to produce detectable 
pulsations, even in sear ches targeted at their known or 
suspected spin periods (I Gotthelf et al.ll2004l ; I Tarn et all 
I2006t iMereghetti et al.ll2006bf ). These objects, and pos- 
sibly the newly-identified magnetar IE 1547.0-5408, rep- 
resent a class of magnetars from which pulsations could 
only be detected intermittently, or perhaps not at all, in 
a search like ours. 

2.3.2. Newly-Identified Sources 

Four sources produced periodic signals that have not 
been previously reported. In Table [1] we have listed ba- 
sic information about each source and signal. Figures [2] 
and [4] contain the Fourier power spectra in which the 
candidate signals were discovered. The power spectra 
provided initial estimates of the pulse periods. We then 
refined the periods by (1) computing pulse profiles from 
non-overlapping intervals between 1000 and 10,000 s long 
(depending upon the number of counts from the source), 
(2) measuring their phases by fitting a sinusoid to the 
data, and (3) modeling the differences between the as- 
sumed and measured phases using first-order polynomi- 
als. The pulse profiles are displayed in Figures [3] and [5[ 
Below we remark on some peculiarities to each source: 
CXOU J174728. 0-321445: The «5000 s period from this 
source is longer than 10% of the exposure time, so it 
would not have been identified in our search were it not 
for the facts that the signal has a strong harmonic near 
2500 s. The source is detected in two observations. In 
both observations, it lies near the edge of a detector. The 
707.4 s dither period is clearly detected in the data, but 
the «1000 s dither period (specifically, 1010.9 s in ob- 
servation 4567 and 989.9 s in observation 4566) is not 
seen because of the orientation of the detector. In addi- 
tion, signals are detected at the difference between the 
frequencies of the dither and the 5000 s signal. The pe- 
riods from the phase connection analysis differ at the 2a 
level, so this signal may not be strictly coherent. 
CXOU J 182531.4-144036: This source was identified for 
its 780 s period. The signal is unrelated to the dither 
periods, which in this observation were at 999.7 s and 
707.6 s. Oddly, a «5000 s period is also present in this 
source. We have confirmed that this 5000 s period does 
not appear in most sources observed with Chandra, and 
the profile of the 5000 s signal is more sinusoidal in this 
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TABLE 1 
Candidate Periodic Signals 



Source 




ObsID 


Detector 


Counts 


-?exp 
(ks) 


Period 
(s) 


-Pmeas 


CXOU J174728.0- 


-321445 


4567 


ACIS-S 


2076 


42.1 


4910±20 


199 






4566 


ACIS-S 


1489 


28.3 


4790±50 


81.9 


CXOU J182531.4- 


-144036 


5341 


ACIS-I 


548 


18.0 


780±3 


36.4 














5000±100 


42.8 



XMMU J124429.7- 


-630407 


010948101 


pn 


491 


49.0 


475.0±0.6 


21.3 








MOS-1 


187 


52.0 


474.6±0.8 


14.4 








MOS-2 


162 


52.0 


474.5±1.0 


3.5 






010948401 


pn 


403 


39.3 


475.1±0.5 


23.5 








MOS-1 


145 


41.7 


476.7±1.7 


11.1 








MOS-2 


141 


41.1 


473.9±1.2 


11.8 


XMMU J185330.7- 


-012815 


0201500301 


pn 


20637 


20.0 


238.2±0.1 


74.2 








MOS-1 


6781 


20.1 


238.4±0.2 


39.1 








MOS-2 


6953 


20.2 


238.4±0.2 


26.3 



Note. — The columns are as follows: the source names; the identifiers of the observations 
in which the signals were found; the detector with which the sources were observed; the total 
number of counts extracted for the source (includes background); the exposure times of the 
observations; the periods, which were computed by tracking the phases of the oscillations; 
and the power which the source was identified in the initial power spectrum. 
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Fig. 2. — Power spectra from sources with periodic variability in 
archival Chandra observations. The left panels display power spec- 
tra from the two observations of CXOU J174728. 0-321445, and 
the right panel displays it from CXOU J182531.4-144036. The 
downward-pointing arrows denote the fundamental and two har- 
monics of the ~700 s and KilOOO s dither periods for the satel- 
lite. The ^700 s dither period, along with beat periods between 
that and the 5000 s signal, are evident in the data from CXOU 
J174728. 0-321445, because the source lay at the edge of the detec- 
tor. 



source than in CXOU J174728.0-321445. However, only 
4 cycles of the 5000 s signal were covered by the obser- 
vation, so this signal may be low-frequency non-periodic 
noise that randomly produced a strong peak in the power 
spectrum. 

XMMU J124429.7-630407: A 475 s period was detected 
from this source in 5 of 6 trial power spectra from two 



Fig. 3. — Profiles of the periodic signals detected in archival 
Chandra observations. Two cycles are repeated in each pane. The 
left panels display the profiles of the ss5000 s signals from the two 
observations of CXOU J174728.0-321445. The right panels display 
the profiles of the 780 s (top panel) and 5000 s (bottom panel) 
signals from CXOU J182531. 4-144036. The period is printed at 
the bottom of each panel. 



different observations. The source is faint, and the non- 
detection with the MOS2 in observation 010948101 can 
be ascribed to unfavorable noise. 

XMMU J 185330.7-012815: This source, also known as 
AX J1853. 3-0128, exhibited a 238 s period in all three 
detectors in the one observation of it. The analysis of 
the phases of this source indicate the signal is not coher- 
ent, as the phases vary by «0.1 cycles with a time scale 
of «5000 s. The periodic signal from this source was 




Fig. 4. — Power spectra from sources with periodic variability in archival XMM-Newton observation. From left to right, we display power 
spectra from the first and second observations of XMMU J124429.7-630407, and from the one observation of XMMU J185330.7-012815. 
From top to bottom, we display data from the EPIC-pn, MOS1, and MOS2. 
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Fig. 5. — Profiles of the periodic signals detected in archival Chandra observations. Two cycles are repeated in each panel, and data from 
all of the pn, MOS1, and MOS2 have been combined. From left to right, we display data from first and second observations of XMMU 
J124429.7-630407, and from the one observation of XMMU J185330.7-012815. 



identified independently by J. Halpcrn and E. Gotthclf, 
who also carried out spectroscopy of its optical counter- 
part, and suggest that it is a cataclysmic variable (private 
communication) . 

We note that it is curious that the period of the sig- 



nal from XMMU J124429. 7-630407 is within 2a of being 
a factor of 2 longer than the period of the signal from 
XMMU J185330.7-012815. We have not been able to 
identify any systematic effect that might explain both 
signals. In our search, we did not find similar significant 
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signals from any of the thousands of other sources that 
we examined. The signals are unlikely to be detector 
artifacts, because they are detected in the pn and both 
MOS cameras. Having failed to identify any causes in- 
trinsic to the spacecraft, on-board computers, or data 
processing, we believe that they are astrophysical, and 
that the factor of two difference between them is simply 
a coincidence. 

3. DISCUSSION 

The periodic signals that we have detected from 
previously-unidentified sources all have periods longer 
than 100 s. Although there is a chance that some of these 
systems are neutron stars, we expect that most of them 
belong to the much-larger population of magnetically- 
accreting white dwarfs, either pola r s or intermediate po- 
lars fe.g.. iNorton fc Watson! 119891: ISchwarz et all 120021 ; 
iRamsav fc Cropper! l2003UMuno et al.ll2003f) . Accreting 
white dwarfs typically have luminosities of 10 30 — 10 32 
erg s -1 , so that given their count rates, they could be 
located as close as 50 pc (for XMMU J185330. 7-012815) 
or as far as 5 kpc (for XMMU J124429. 7-630407). As 
members of the local Galactic neighborhood they will be 
distributed over the entire sky, not just within |6|<5°, 
our targeted survey of the Galactic plane is not particu- 
larly efficient at finding cataclysmic variables. Therefore, 
we do not discuss them further. 

In the context of a search for magnetars, it is notable 
that no new source was found to exhibit periodic vari- 
ability with periods in the range of known magnetars 
that are pulsed in X-rays, 5-12 s. To understand the 
lack of detections of obvious candidate neutron stars, 
we need to compute the fraction of the Galaxy that 
was covered by our survey. If the properties of mag- 
netars as a population were bettern known, we would do 
so by assuming distributions for their luminosities and 
pulse amplitudes, and carry out a maximum-likelihood or 
Monte-Ca rlo calculation to model the observed popula- 
tion (e.g., iFaucher-Giguere fc Kaspill2006t lLorimer et al.1 
2006). However, the intrinsic distributions for the lu- 
minosities and pulse fractions of magnetars are poorly 
constrainted. Some guidance can be obtained from mod- 
els that explain the pulsations as originating from single 
bright spots on their surfaces. lOzell (2002) find that the 
fractional amplitudes of pulsations are largest when the 
hot spot is located on the equator, and viewed from the 
equator. However, the amplitude drops by «50% when 
the spot and viewer have a latitude >50°. Therefore, we 
roughly estimate that any given magnetar will be easily 
detectable over «65% of the sky. 

Unfortuntely, as we describe below, the luminosities 
of magnetars are highly variable, and cannot be pre- 
dicted by first principles because the mechanism causing 
the variability is n ot understood (e.g.. lWoods et al.ll2004h 
iMuno et al.ll2007T ). Therefore, in the following we only 
present some representative cases. The total Galactic 
populations for our fiducial examples are then calculated 
in two steps, first computing the depth along the line-of- 
sight through the Galaxy to which our observations were 
sensitive, and second estimating the fraction of the stel- 
lar mass in the Galactic spiral arms that was enclosed by 
our observations. 

3.1. Depth of the Survey 



We compute the depth (D) through the Galaxy that 
each observation was sensitive for any given luminosity 
(Lx) and limiting total number of counts (Ci; m ) from: 

C Um /T exp = Ua(N s [1, b, D])^^, (4) 

where £,det(Nn[l,b, D]) is the conversion factor between 
flux and count rate for each detector. This factor de- 
pends additionally upon the Galactic absorption Ah, 
which in turn is a function of the distance D to a source 
and its position in the Galactic plane, b. We com- 
puted /(Ah) for both the Chandra ACIS-I and for the 
XMM-Newton EPIC pn behind the medium filter using 
the Portable, Interactive Multi-Mission Simulator, 9 as- 
suming several trial spectra and a range of Ah- We have 
neglected some other factors that do not significantly af- 
fect our results, the choice of filters for XMM-Newton 
(an «5% effect on £), and hydrocarbon build-up on the 
ACIS (negligible for sources with Ah>10 22 cm -2 ). 

Vignetting as a function of offset from the aim point is 
a small effect for Chandra, reducing the count rate from 
a source by ~20% at an offset of 8'. We have neglected 
vignetting in our Chandra observations. However, it is 
a large effect for XMM-Newton, reducing the count rate 
by 50% at an offset of 10'. Therefore, we have accounted 
for vignetting in our XMM-Newton observations by re- 
ducing the flux-to-counts conversion £ by 33% from the 
on- axis value obtained from PIMMS, which is the mean 
reduction over the inner 10'. We estimate that our sim- 
ple treatment of the vignetting introduces an uncertainty 
of <15% on the depth of the survey. 

The spectra of magnetars are generally described as 
the sum of a blackbody component with a temperature 
of A;Tss0.6 keV and a soft power-law tail with photon 
index T«3 (e.g.. I Woods fc Thompsonll2006f K The overall 
spectrum can be roughly approximated as a T=2 power 
law, so we take that as our fiducial spectrum. Choosing 
instead a softer T=3 power law results in values of D that 
are up to 30% smaller (depending upon the absorption 
through the line of sight; see below), while choosing a 
harder A:T=0.6 keV blackbody increases D by up to 20%. 

We estimated the absorption as a function of distance 
along different lines of sight using models for Galactic 
optical and infrared extinction. For most of the Galac- 
tic plane, we linearly interp olated visual ex t inctio n val- 
ues (Ay) from the model of iDrimmel et al] ((2003), and 
converted them to K b and extinction values u s ing the 
relati on A K /A V = 0.11 (jRieke fe Lebofskvil 19851: iMathisI 
1990). However, withi n the central 25° of the Galaxy the 
model extinction from IDrimmel et alJ (l2003f) was signif- 
icant ly lower than that observed fe.g.. lLaunhardt et al.l 
2002), so instead we interpo lated values of the K s band 
extinction from the table in iMarshall et alJ (|2006l ) . We 
then converted the K band extinction (ignoring the 
5% difference in Ak and Ak s ) int o a column den- 
sity using A H = 1.6 x 10 23 Aj f cm 2 (jRieke fc Lebofskvl 
[T98llPredehl fc Schmittl [l995). Absorption through the 
Galactic plane can reduce the flux observed from a source 
by up to a factor of 30 compared to the value without 
absorption. We tested different values for N-r/Ak in the 
ran ge (1.2 : — 2.1) x 10 22 cm 2 (repre senting values from, 
e.g. "lGlassl[l999l : lTan fc Draine!l2004D . and found that the 

9 http: //heasarc .gsf c .nasa. gov/docs/sof tware/tools/pimms .html 



A Search for Magnetars 



9 



uncertainty in D introduced by our choice of absorption 
model is «15%. 

The depth of our survey depends strongly upon the as- 
sumed luminosity and the fractional root-mean-squared 
(rms) amplitude of the pulsations (A ms ) for which we 
are searching. The rms amplitude enters consideration 
because it determines Cn m : 



100 r 



1.13 



2P • 

Zj± Slg 



lim 



1/2 



7T V 

2 



^Nyq 



(5) 



where P s j g is the intrinsic power of a detectable signal, 
the factor of 1.13 is an average correction that accounts 
for the fact that signals will often fall between indepen- 
dent Fourier bins, and the sine term takes into account 
the attenuation of the power of a signal as its frequency 
y ap proaches the Nyquist value j^Nyq (Vaughan et al. 
1994). The intrinsic power P s j g detectable in a search 
must be determine d using the distr i bution of noise pow- 
ers, as described in IVaughan et all |l994). For a search 
threshold of P me as=32.8, the expected value (50% con- 
fidence) is P s i g =31.8, and the power detectable in 90% 
of trials is P s j g =48.0. In the following we will use the 
90% confidence level, P s i g =48.0. The value of D com- 
puted for P s i g =31.8 is wlO-20% larger than the value for 
P sig =48.0. 

For our depth calculation, we use two extremes as ex- 
amples. The first is an easily-detectable magnetar. As 
our model for this, we use a faint pulsar with a high 
pulse frac tion, such as CXOU J164710.2-455216 in West- 
erlund 1 (jMuno et alJ 12006). which had a luminosity of 
L x =3 x 10 33 erg s _1 (0.5-10 keV) and a fully-modulated 
sinusoidal pulse profile with A rms =0.71. Pulsations like 
this could be detected in 90% of trials with as few as 
wl20 photons. The depth probed by any given ob- 

1/2 

servation scales as L x and linearly with A lms , so the 
easily-detectable example is equivalent to a bright mag- 
netar like SGR 1900+14, with L x =10 35 erg s" 1 (0.5- 
10 keV) and a sin usoidal pulse profile with A rms =0.12 
([Hurley et al. I ll 999). In the top panel of Figured! we plot 
the depth at which each observation would have been 
sensitive as a function of exposure time for our easily- 
detectable example. Only 10% of the Chandra observa- 
tions and 20% of the XMM-Newton observations probe 
the entire depth of the Galactic plane to a distance of 24 
kpc from Earth. 

The other extreme is a barely-detectable pulsar, with 
a low luminosity Lx=3 x 10 33 erg s _1 and a small pulse 
fraction A Ims =0.15. This example is chosen to have been 
just detectable in observations comparable to those of 
XTE J1810-197 in quiescence, for which L x =2 x 10 33 erg 
s _1 (0.5-10 keV), and from whic h pulsations were not de- 
tected to a limit of A rms <17% (jGotthelf et al.ll2004h . A 
source with A rms =0.15 would require 2800 photons to be 
identified in 90% of trials. Consequently, a search sensi- 
tive to our barely-detectable pulsar would cover a factor 
of «5 less depth than a search for easily-detectable mag- 
netars like the one in Wester lund 1. We display the depth 
of a survey for a barely-detectable magnetar in the bot- 
tom panel of Figure [6] None of the archival observations 
would probe the entire Galaxy when searching for our 
barely-detectable pulsar. 

3.2. Fraction of Young Stars Enclosed by the Survey 




Open symbols: 
A r „„ = 0.15 



Exposure (s) 

Fig. 6. — The depth through the Galaxy that each observation 
sampled, plotted as a function of its exposure time. Triangles illus- 
trate Chandra ACIS observations, and circles XMM-Newton EPIC- 
pn observations. The EPIC-pn has a larger effective area than 
ACIS, so the XMM-Newton observations are systematically deeper 
for a given exposure. Filled symbols indicate the depth at which 
we are sensitive to an easily-detectable pulsar with Lx = 3 X 10 33 
erg s — 1 and j4 rma =70%. Open symbols indicate the depth at which 
we are sensitive to a barely-detectable pulsar with Lx = 3 X 10 33 
erg s and ^4 rms =15%. The dashed line at a depth of 24 kpc de- 
marks observations that proble the entire Galaxy; larger depths are 
no longer meaningful because they extend beyond the distribution 
of young stars. 



In order to estimate the number of magnetars in the 
Galaxy and their birth rates, we need to know the frac- 
tion of the Galaxy that we have meaningfully surveyed 
for magnetars. Magnetars have massive progenitors: 
neutron stars form from stars within inital masses >8Af Q 
(jHeger et al.l 12003). and there is evidence that mag- 
netars form from >30 stars (|Gaensler et al.l 120051 : 
iFiger et all 120051: iMuno et all |2006l) Th e lifetimes of 
magnetars are thought to b e ~10 4 yr (|Kouveliotou et all 
1999: iGaerlsler et al.l fl999h . so even if some magnetars 
receive kicks of ^1000 km s -1 , they will be found within 
~10 pc of their birth place. Most massive stars are con- 
centrated in the Galaxy's spiral arms, so we have calcu- 
lated the fraction of the stellar mass in the spiral arms 
that our observations have covered. 

We use the model for the st ellar mass in the spiral arms 
from IWainscoat et alJ |l992).The locations of the arms 
are defined as a logarithmic spiral: 



9(R) = alog(R/R mia ) + e is 



(6) 



where R is the radial distance from the Galactic center, 
9 is the azimuthal angle in radians (with 6=0 along the 
line connecting the Earth to the Galactic center, and 
O<0<7r for sinZ<0), Pt m ; n is the radial distance at which 
the arms start, m i n is the angle at which the arms start, 
and a is the winding constant. We define four main 
spiral arms, and one smaller, "local arm." Each arm 
is assumed to extend through an angle # cx t, out to a 
maximum radial distance of P max =15 kpc. The values 
of a, -Rmin, fl mim and flext are listed in Table O taken 
directly from IWainscoat et al.l ([1992). The parameters 
assume the Earth is 8.5 kpc from the Galactic center. 

The stellar density along the spiral arms is given by an 
exponential distribution with radius R and height above 
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TABLE 2 
Spiral Arm Parameters 



Arm 


a 


-f^min 




#ext 




(radians) 


(kpc) 


(radians) 


(radians) 


1 


4.25 


3.48 


0.000 


6.00 


1' 


4.25 


3.48 


3.141 


6.00 


2 


4.89 


4.90 


2.525 


5.47 


2' 


4.89 


4.90 


5.666 


5.47 


L 


4.57 


8.10 


5.847 


0.55 



the plane z (jWainscoat et al.l [l992) : 

p = Po cxp(-R/h-\z\/h z ), (7) 

where /i=3.5 kpc is the radial scale length, h z —90 pc 
is the scale height of the youngest stars in the Galactic 
plane, and po is a normalization for which the exact value 
is unimportant, because we will divide by the total mass 
to obtain the fraction of the arms encompassed by the 
survey. The density of the spiral arms is assumed to be 
constant for any given R within ±375 pc of the center 
of the arm, and zero outside that radial extent. This 
model for the spiral arms is illustrated schematically in 
Figure [7] In that figure, we also plot the depth through 
the Galaxy for which we would detect easily- and barely- 
detectable pulsars in 90% of trials (P s i g =48.0) for each 
of our archival observations. The 50% completeness level 
would incorporate about 20% more of the Galaxy, but 
using it would add the complication that a larger fraction 
of magnetars with Lx and A rms above our stated limits 
would be missed. 

For each observation, we integrated the mass enclosed 
along the line of sight out to the depths defined in §3.1. 
We assumed a field of view of 17' x 17' for Chandra, which 
corresponds to the full ACIS-I array. We assumed a field- 
of-view of 20' x 20' for XMM-Newton, which corresponds 
to the region where vignetting still provides >50% of the 
on-axis count rate for a source. The fractional mass of 
the spiral arms enclosed by each observation is plotted 
as a function of Galactic longitude in Figure for the 
cases of easily- and barely-detectable pulsars. 

To compute the total fractional mass enclosed by the 
survey, we identified duplicate observations as those 
within 12' of each other, and kept only the deepest of 
the duplicates. Then, we summed the fractional masses 
of the unique observations in Figure[5] We estimated un- 
certainties on the total fractional masses by comparing 
the results of the calculations made with different values 
for the fiducial spectrum and for the conversion between 
infrared extinction and X-ray absorption (§3.1). We find 
that our survey is 90% complete (P s i g =48.0) for easily- 
detectable pulsars (Lx=3 x 10 33 erg s _1 , A rms =71%) for 
5% of the Galactic spiral arms under our standard model. 
The XMM-Newton and Chandra observations are about 
equally efficient, with Chandra surveying 2% of the spi- 
ral arms, and XMM-Newton 3% of the spiral arms. Us- 
ing reasonable alternative spectra and absorption values 
described in §3.2, the total mass fraction surveyed can 
range between 3% and 7%. For barely-detectable pulsars 
(Xx=3 x 10 33 erg s _1 , A rms =15%) our survey is complete 
for 0.4% of the spiral arms, with a range of 0.2-0.6% if we 
choose different input values. XMM-Newton and Chan- 
dra each surveyed 0.2% of the spiral arms on their own. 




X (kpc) 

Fig. 7. — The model for the Galactic spiral arms, viewed from 
above the north Galactic Pole, overlaid by the depth of the ob- 
servations in th e surv ey. The black lines are the model from 
Wainscoat et al. (1992). The blue lines represent Chandra obser- 
vations, and the red XMM-Newton observations. The top panel 
illustrate the depth to which we are sensitive to easily-detectable 
pulsars, and the bottom to barely-detectable ones. 



We note that this model for the spiral arms does not 
include any young stars in the cen tral 150 pc of the 
Galaxy, which has been surveyed by IWang et all (|2002f ) 
and Muno et al. (in prep). For the easily-detectable mag- 
netar case, our archival survey of this region reaches the 
Galactic center (Fi gure [71). This region co ntains ~1% of 
the Galactic mass ()Launhardt et al]|2002f). but the star 
formation rate is still under debate. iFiger et all (|2004D 
modeled the stellar population observed in the infrared 
there, and concluded that the star formation rate is ~1% 
of the Galactic value. However, it is pos sible that sta r 
formation is skewed toward massive stars (|Morrisl [l993). 
This is suggested by indirect measurements of the Lya 
flux in the region, which could be 10% of the total Galac- 
tic value (|Cox fc Laureiid 119891 : IFiger et al.lll999t ). De- 
pending upon the true star formation rate, the central 
2°xl°of the Galaxy could contribute more than half of 
the massive stars studied in this project. 

3.3. Constraints on the Population of Magnetars 

The main uncertainty in the birth rate of magnetars 
in the Galaxy is what the "typical" magnetar looks like. 
Most known examples are easily-identifiable in Chandra 
and XMM-Newton observations, by which we mean that 



A Search for Magnetars 



11 



o- 3 

cr 4 

o- 5 

o- 6 

cr 7 
o- 8 



iChandra 
»XMM 




A r „« = 0.71 



f. 



-100 





(deg) 



100 




= 0.15 



(deg) 



Fig. 8. — The fraction of the mass of the spiral arms surveyed 
by each observation, plotted as a function of Galactic longitude. 
Triangles illustrate Chandra ACIS observations, and circles XMM- 
Newton EPIC-pn observations. The EPIC-pn has a larger effective 
area than ACIS, so the XMM-Newton observations are system- 
atically deeper for a given exposure. The top panel, with filled 
symbols, indicates the mass surveyed for which we are sensitive 
to an easily-detectable pulsar, with Lx = 3 X 10 33 erg s _1 and 
A rms =70%. The bottom panel, with open symbols, indicates the 
mass surveyed for which we are sensitive to a barely-detectable 
pulsar with Lx = 3 X 10 33 erg s _1 and A rms =15%. 



they are either: (1) luminous (Lx^lO 35 erg s x ) with a 



modest pulse fraction (A rms >13%) or (2) faint (Lx~3 x 
10 33 erg s -1 ) with a large pulse fraction ( J 4 rms «70%). 
Indeed 7 of 12 Galactic examples were identified 
blindly in our search (§2.4). However, pulsations 
are only intermittently-detect able from three transien t 
magnetars — SG R 1627-41 dMereehetti et alJ [2006b]). 
AX J1845.0-0258 (iTam et all 120061). and XTE J1810- 
197 (|Gotthelf et ; all 120041 ). A fourth, IE 1547.0-5408, 
has recently b een identified as a magnetar through its 
radio emission (|Camilo et al.ll2007D. but X-ray pulsation s 
have not yet been confirmed (|Gelfand fc Gaenslerll2007l) . 
The limits on their luminosities and pulse fraction serve 
as guides for our barely-detectable pulsars. Therefore, 
there could be a significant number of magnetars that 
can only be identified in a blind search intermittently, 
when they produce outbursts. 

For the purposes of this discussion, we will divide 
magnetars into three types: (1) the standard AXP- 
like sources, which are persistently bright and can be 
detected with pointed observations at any time, (2) 
the SGR-like magnetars, which can be identified by 
wide-field gamma-ray burst monitors from anywhere in 



the galaxy whenever they produce outbursts consisting 
of multiple soft gamma-ray bursts or occasional giant 
gamma-ray flares, and (3) the transient AXPs, which 
can only be identified by pointed observations when 
they produce outbursts reaching a luminosity of >10 35 
erg s _1 (0.5-10 keV). This is not to suggest that these 
groups are mutually exclusive. For instance, the stan- 
dard AXP IE 2259+586 exhibits both short time scale 
bursts that appear similar to (although much fainter 
than) SGR bursts, and variability in its mean inten- 
sity that could be interpreted as a transient outburs t 
with with time scales of months ([Woods et all [2004) . 
Moreover, in addition to its soft gamma-ray bursts, SGR 
1627-41 also exhibits large (factor of 20) va riations in its 
persistent X-ray luminosity between bursts (IWoods et al.l 
ll999UHurlev et al.ll2000HMereghetti et al.ll2006al ). which 
could be interpre ted as an outburst l i ke those seen from 
trans i ent AXPs (IWoods et all 120011 : iKouveliotou et all 
120031: IWoods et alJ 120061 ). Instead, we use these clas- 
sifications to highlight the relative ease or difficulty with 
which different magnetars could be detected, and thereby 
to answer the question, If there are many more magne- 
tars in the Galaxy, what must they look like? 

In the subsections below, we use the results of sur- 
veys most sensitive to each class of object. For standard 
AXPs, we use our own survey of Chandra and XMM- 
Newton data. For SGRs, the best constraints are pro- 
vided by all-sky gamma-ray burst monitors. For tran- 
sient AXPs, similar results are found using our own 
archival survey, and past surveys with RXTE and ASCA. 

3.3.1. The Number of Standard AXPs 

The standard AXPs are the easiest population to con- 
strain, because they are easily-detectable throughout the 
Galaxy. Our survey of archival Chandra and XMM- 
Newton observations is 90% complete for finding stan- 
dard AXPs for «5% of the young stellar population in 
the Galaxy. 

Before estimating a total number of standard AXPs, 
however, we need to establish how many known exam- 
ples lie within our random survey of the Galaxy. Only 
one source was identified in an observation in which it 
was not t he target: CXOU J164710.2-455216 in West- 
erlund 1 (Mu no et alJ 120061 ). The other sources de- 
tected in our survey were observed as the targets of 
Chandra and XMM-Newton observations because they 
were previously known, and some thought must be taken 
before they can be considered as part of a random 
sample. Three of the standard AXPS originally were 
identified serendipitously during observations of other 
sources: IE 1048.1-5 937 during Einstein observation s 
of the Carina nebula (|Seward. Charles, fe Smale] [l986). 
IE 2259+586 during Einstein observatio ns of a super- 
nova remnant (|Fahlman fc Gregory|[l98ll) . and IE 1841- 
045 during ASCA observat ions of a supernova remnant 
(jVashisht fc Gotthelflll997l) . In the first case, the sepa- 
ration on the sky between the magnetar and the central 
position of the original target were >0.5°, which makes 
it unlikely that the magnetar would have serendipitously 
fallen in the field of view of a Chandra or XMM-Newton 
observation. In the latter two cases, the magnetars are 
close enough to the supernova remnants that one could 
have discovered them with Chandra or XMM-Newton ob- 
servations, if one makes the reasonable assumption that 
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3.3.2. The Number of SGRs 
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Fig. 9. — The number of SGRs that could be present in the 
Galaxy, as a function of the rate at which year-long outbursts (con- 
sisting of multiple, <1 s soft gamma-ray bursts) recur. The solid 
lines shows the most likely number, given that 3 Galactic SGRs 
have been identified in the last 25 years, and the dashed lines show 
the 90% confidence limit on that number. 



the supernova remnant would be a sufficiently compelling 
target to have been observed. Therefore, we claim that 3 
known magnetars lie in the random sample of the Galaxy 
that we surveyed. 

Knowing that there are 3 easily-detectable magnetars 
that can be found in our search of 5% of the young stars 
in the Galaxy, we can determine the most likely total 
number of such objects in the Galaxy using the binomial 
distribution. The probability that n pulsars out of a total 
population N would lie in a fraction / of the Galaxy is: 



prob(n|/,iV) = /»(l-/)( Ar -") 



Nl 



nl(N-n)\ 



(8) 



This can be inverted using Bayes' theorem to give the 
total number of magnetars N, given that n—3 magnetars 
are found with the fraction /=0.05 of the Galaxy: 



prob(iV|/, n) 



prob(n|/, N) 
ESUprobM/.W)' 



0) 



The most likely value and 90% confidence interval for 
the total number of easily-detectable magnetars is 591^2- 
Given that the number of known Galactic magnetars is 
12, there could be 47+32 waiting to be discovered. 

The Galactic rate of core-collapse super- 
novae is ~0.01-0.04 yr" 1 dCappallaro et ail 119931 : 
Ivan den Bergh fc McClurel Il994f) . and the radio pulsar 
birth rate has been recently estima ted to be between 
^0.014 yr" 1 (lLorimer et all I2006D and «0.03 yr" 1 



(jFaucher-Giguere fc Kaspil I2006D. If magnetars have 
typical lifetimes of ~10 4 vr (|Kouveliotou et alj 119991 : 
iGaensler et al.lll999l ). then with somewhere between 27 
and 160 in the Galaxy, the magnetar birth rate would 
be between 0.003 and 0.016 yr _1 . The lower bound 
is considerably smaller than the pulsar birth rate, as 
would be expected if magnetar s descend from >30M© 
stars (e.g.. IGaensler et a.l.ll2005h . However, at the upper 
range, magnetars that resemble the standard AXPs 
could have a birth rate nearly equal to that of radio 
pulsars. 



An SGR can be identified if it produces multiple 
bright soft gamma -ray bursts (>10 40 erg s -1 , e.g., 
lAptekar et al.ll200l[ h Over the course of «25 years of 
monitoring, 3 Gal actic SGRs have been confi rmed (and 
one in the LMC I Woods fc Thompson! l2006h. 1 Their 
bursts typically last ~0.1 s (Go giis et al.l 120011) . so they 
are bright enough to be detected from anywhere in the 
Galaxy. Therefore, to estimate their total number, we 
only need to determine the rate at which SGRs produce 
detectable bursts. 

For any monitoring period of length T, if we define an 
"outburst" as containing multiple individual soft gamma- 
ray bursts, an SGR will be identified if it produces one 
or more outbursts. According to Poisson statistics, if the 
outburst rate r (where we take r to have units such that 
the expected number of outbursts in a time T is rT), the 
chance that any given source is detected will be 

prob dct = 1 - er rT . (10) 

Given an ensemble of n sources, the probability that m 
sources will be detected (and n — m not detected) is 

prob(m|n,rT) = (1 - e -»Tjm( e -rrj(n-m) 



m\(n — m)\ 

This is simply a binomial distribution, which can be in- 
verted using Bayes' theorem to obtain prob(n|m, rT) in 
the same way as Equation [9l 

We can receive some guidance about the rate at 
which SGRs are active by noting that, out of 25 years 
of monitoring, SGR 1806-20 has been active for sa7 
years, SGR 1900+14 for w4 years, SGR 0526-66 (in the 
LMC) for ^2 years, and SGR 1627-41 for »2 months 
(|Aptekar et ali 120011: IWoods fc Thompson! I2006D . The 
bursts tend to be associated with periods of activity last- 
ing months to years. It is clear that the activity level is 
variable from source to source, so we consider a range of 
rates. 

Given m—3 and T=25 years, in Figure [9] we plot for 
a range of r the most-likely value and 90% confidence 
limit for n. If on average one outburst is expected ev- 
ery 10 years, which is near the mean of the sample that 
has already been identified, then with 90% confidence we 
should have already detected all of the Galactic SGRs. 
However, if an SGR only becomes active once per 100 
years, there could be I3I7 8 m total in the Galaxy. Our 
best estimate is therefore that the number of Galactic 
SGRs is much smaller than that of standard AXPs, un- 
less SGRs are active less frequently than once per several 
hundred years. 

3.3.3. The Number of Transient AXPs 

Transient AXPs would be the most difficult to iden- 
tify, and could comprise the majority of the magnetars. 
We define an AXP as transient if it has a luminosity of 
<10 33 erg s _1 (0.5-10 keV) for most of its current life, 
but increases in luminosi ty to ^10 35 erg s _1 for durations 



of on order a year (e.g., iGotthelf et all 12004 iTam et all 



10 In addition, one candidate SGR has been pro posed based on 
the identification of two bursts (Clinc et al. 2000), and hundreds 
of individual SGR-like bursts have been identified all over the sky 
(see www.srl.berkeley.edu/iph3). 
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Fig. 10. — The number of transient AXPs that could be present 
in the Galaxy, as a function of the rate at which year-long outbursts 
(in which the X-ray luminosity is at least 10 35 erg s — 1 ) recur. The 
number of Galactic examples is calculated based on observations 
with ASCA, RXTE, Chandra, and XMM-Newton, which surveyed 
~2% of the Galaxy each year for 11 years, and discovered two 
transient AXPs. The solid lines shows the most likely number, 
and the dashed lines show the 90% confidence limit. 



2006). We define the "outburst" as the year-long period 
when the AXP is luminous. Three transient AXPs are 
known (XTE J1810-197, CXOU J164710. 2-455216, and 
IE 1547.0-5408), along with one candidate (AX J1845.0- 
0258). 

Our archival Chandra and XMM-Newton survey pro- 
vides the most sensitive search for transient AXPs in 
quiescence. Nonetheless, when looking for a pulsar with 
L x =3 x 10 33 erg s" 1 (0.5-10 keV) and A rms =15%, our 
survey is only 90% complete for «0.4% of the Galac- 
tic young stellar population. The simple calculation 
(Eq. [9|) assuming no new barely-detectable magnetars in 
/=0.005 of the Galaxy, 11 implies that with 90% confi- 
dence there could be up to 540 waiting to be discovered. 
Unfortunately, this estimate is not a strict upper bound, 
because quiescent AXPs could be less luminous or have 
lower pulsed fractions than we have assumed. 

An AXP in outburst could be detected over a larger 
fraction of the Galaxy, but it still requires an observation 
with a pointed X-ray observatory. To constrain their 
numbers based on AXPs in outburst, we must account 
for both the fraction of the young stars in the Galaxy 
that can be surveyed by a given observatory (Eq. [8|), 
and the chance that AXPs in the region will produce an 
outburst (Eq [TTj) . We assume that there is a population 
of N magnetars in the Galaxy with an outburst rate r, 
and that we searched for them with a survey of duration 
T that covered a fraction / of the magnetar birth places. 
Then, the probability that m magnetars will be found is 
the joint probability that n out of N magnetars will lie 
in the survey region, and that m out of n will produce 
outbursts: 

prob(m|rT, /, N) = 

11 CXOU J164710. 2-455216 is transient, but was already counter 
among the easily-detectable examples, because it has a fully- 
modulated pulse profile that is not characteristic of the other tran- 
sient AXPs. 



x/"(l - 



JV! 



n\(N -n)V 



(12) 



This can be inverted using Bayes' theorem, in order to 
derive the probability that there are N magnetars given 
that m are found with outbursts at a rate r in a frac- 
tion / of the Galaxy, prob(7V|m, rT, /), as in Equation[9l 
Here, because we are dealing with transient sources, / 
represents the fraction of the young Galactic population 
searched per year, which is the duration of a typical out- 
burst. 

We can estimate the rate of outbursts from tran- 
sient AXPs from the four known examples (here we 
assume that the candidate is a genuine example). 
AX J1845. 0-0258 and XTE J1810-197 have each ex- 
hibited only one outburst each, in 1993 and 2002 re- 
spe ctively, which in both cases led to their discover- 
ies llTorii et alj|1998t llbrahim et al]|2004t iGotthelf et all 
I2004t iTam et all 12006ft . In contrast. CXOU J164710.2- 
455216 entered into outburst only a year after it s dis- 
covery (jMuno et aUl2006L 120071 : llsrael et al.ll2007h . and 
IE 1547.0-5408 has brightened to ~10 35 erg s" 1 (0.5- 
10 keV) on at least two occasions in the last 10 years 
(lOelfand & Oaenslerl [20071 : iCamilo et al1l200l . Taken 
together, these results suggest that the mean recurrence 
time is about once every 8 years. As we did for SGRs, 
we consider below a range of rates around this value. 

To compute the fraction of the young Galactic popula- 
tion that has been searched for transient AXPs, we need 
to take into account how the known examples were dis- 
covered. Two transient AXPs were discovered serendip- 
itously: XTE J1810-197 in an R XTE observation of 
SGR 1806-20 (jlbrahim et al.ll200l . and A X J1845.0- 
0250 in an ASCA observation of Kes 75 (|Torii et all 
1998). We consider these to be useful examples for es- 
timating the total Galactic population. The other two 
probably would not have been idcnitified as transient 
AXPs were it not for unusual circumstances. The out- 
burst from CXOU J164710.2-455216 was identified be- 
cause it produced a hard X-ray burst that was detected 
by Swift, but that burst probably would have been ig- 
nored if the source had not b een previously identified 
in quiescence ([Israel et al.ll2007T) . We have taken into ac- 
count objects like CXOU J164710.2-455216to the Galac- 
tic population in §3.3.1, so we do not include it here. 
IE 1547.0-5408 was positively dentified as a magnetar 
through radio observations of a variable, point-like X-ray 
source i n a supernova remnant that was studied delib - 
erately ijGelfand fc Oaensled 120071: ICamilo et al]l200l . 
X-ray pulsations have not yet been confirmed from this 
source, so we do not consider it here. 

Based on the above considerations, we now calculate 
the fraction of young stars that RXTE, ASCA, Chan- 
dra, and XMM-Newton surveyed per year during their 
lifetimes, in the same manner as in §3.1-3.2. We as- 
sume that the target magnetar would have Lx=10 35 
erg s _1 (0.5-10 keV), and a net number of counts de- 
termined by the observations in which XTE J1810- 
197 and AX J1845. 0-0258 were identified. For RXTE 
XTE J1810-197 was detected with «6 counts s _1 in a 
2.6 ks observation with the Proportional Counter Array 
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(PCA; llbrahim et aD 12004). so we assume that a pul- 
sar would require 1.6 x 10 4 total counts to be identi- 
fiable. We take the PCA field-of-view to be 50'x50', 
and use a count-to-flux conversion that is 33% of the on- 
axis value. We find that the RXTE PCA has surveyed 
«2% of the Galaxy down to this sensitivity during the 
last 11 years, and ~0.7% in any given year. For ASCA 
AX J1845.0-0258 was detected wi th ftilOOO net cou nts in 
each Gas Imaging Spectrometer (|Torii et al.lll998T) . Vi- 
gnetting drastically reduces the count rates from sources 
off-axis, so we assume a field-of-view of only 20'x20', and 
use a flux-to-count conversion that is 33% smaller than 
the on-axis value. We find that ASCA observations were 
sensitive to bright magnetars over «3% of the Galaxy 
during its 7 year lifetime, and ~1% of the Galaxy per 
year. Finally, we find that Chandra surveyed «0.5% of 
the Galaxy per year over the last 7 years, and that XMM- 
Newton surveyed «1% of the Galaxy per year over the 
last 6 years. Taken together, and accounting for over- 
laps, we find that over the last 11 years (since the launch 
of RXTE), «2% of the Galaxy has been surveyed by a 
pointed X-ray observatory each year. 

In Figure ITOl we plot the constraints on the total num- 
ber of transient AXPs, given that two were identified 
during observations of /=0.02 of the young Galactic pop- 
ulation each year for T=ll yr. Transient AXPs are prob- 
ably the largest undiscovered population of magnetars. 
For instance, if the average transient AXP produces out- 
bursts at a rate of one per 10 years, then there could 
be 190T^io °f them in the Galaxy. The upper bound is 
comparable to the number of quiescent magnetars that 
we estimate could be present based on archival Chandra 
and XMM- Newton observations. For a lifetime of 10 4 yr 
(|Kouveliotou et all 119991 iGaensler et al.lll999h . this im- 
plies a brigh rate of 0.008-0.06 yr -1 . The lower bound is 
«50% of the radio pulsar birth rate, whereas the upper 
boun d exeeds the higher estima t e of th e pulsar birth rate 
from IFaucher-Giguere fc Kaspil (|2006l ). The higher rate 
would be large compared to the estimated Galactic su- 
pernova rate. However, it is possible that transient mag- 
netars have significantly longer lifetimes than have been 
estimated base d on the energy budgets of the persistent 
exam ples (e.g., iKouveliotou et all 119991 : IGaensler et alJ 
1999), in which case the birth rates would be correspond- 
ingly lower. 

4. CONCLUSIONS 

We have serached archival Chandra and XMM-Newton 
observations for X-ray pulsars, in order to constrain the 
Galactic population of magnetars. Although we found 
four objects with periodic variability on time scales of 
200-5000 s (Figs. [M]and Table H]), we found no sources 
that were obviously magnetars. The archival obser- 
vations that we used covered a moderate fraction of 
the young stellar populations from which magnetars de- 
scend. Our search was sensitive to the standard, easily- 
detectable AXP-like magnetar (i x =10 35 erg s -1 [0.5-10 
keV] and a sinusoidal pulse profile with A rms =0.12, or 
equivalently Lx=3 x 10 33 erg s _1 and A rms =0.7) with 
90% completeness for ss5% of the mass in the Galactic 
spiral arms. Searching for transient AXPs in quiescence 



Aptekar, R. L., Fredericks, D. D., Golcnctskii, S. V., Il'i 

N., Mazets, E. P., Pal'shin, V. D., Butterwort, P. S., & Cline, T. 
2001, ApJ, 137, 227 



x 10 33 erg s _1 and A rms =0.15), however, we were 
only sensitive over 0.4% of the Galactic spiral arms. 

Based on the number of known Galactic magnetars, 
we then placed constraints on the total population. Our 
archival search placed strict constraints on the number 
of standard, persistent AXPs, which number 59^2 with 
90% confidence. We also placed constraints on the num- 
ber of SGR-like magnetars as a function of the rate at 
which they produce outbursts (Fig. [9]), and find that 
for a likely recurrence rate of 0.01 yr" 1 , there would 
be 13^7 8 in the Galaxy. However, these populations 
may be dwarfed by the number of transient AXPs sim - 
ilar to XTE J1810-197 (see also llbrahim et all l200l . 
No new examples were detected in archival Chandra and 
XMM-Newton observations, so up to 540 quiescent AXPs 
could be present in the Galaxy. We also considered the 
fact that two transient AXPs have been discovered in 
outburst by ASCA and RXTE. The fact that ASCA, 
RXTE, Chandra, and XMM-Newton together have sur- 
veyed 2% of the Galaxy each year for the last 11 years 
allows us to constrain the number of transients as a func- 
tion of their recurrence time (Fig. [TO]) . For a recur- 
rence time of once every 10 years, this calculation sug- 
gests there are between 80 and 580 transient AXPs, with 
90% confidence. Assuming a lifetime of 10 4 yr for persis- 
tent magnetars (|Kouveliotou et ai1ll999t IGaensler et al.l 
1999), their birth rate is at least 10% of that of ra- 
dio pulsa rs, and could equal the radio pulsar birth 
rate (e.g., lLorimer et alj|2006t IFaucher-Giguere fc Kaspil 
2006). The birth rate of transient magnetars could ex- 
ceed those of radio pulsars, unless their lifetimes are sig- 
nificantly longer than 10 4 years. 

The main uncertainty in the number of Galactic mag- 
netars is introduced by the unknown activity level of 
transient examples. Further monitoring of known tran- 
sient AXPs would constrain their duty cycles and re- 
duce this uncertainty, so long as the known examples can 
be taken as representative of all magnetars. Ultimately, 
though, a wide-field monitoring campaign is needed to 
identify more transient AXPs in outburst. This can be 
accomplished in principle in the X-ray band, although 
the sensitivity of current and planned wid e-field moni- 
tors (»2 x 10~ n erg cm" 2 s" 1 , 0.5-10 keV;jLevine et al.l 
ll996tlRemillard et alJl2000b iBraga fc Meji'aH2006h is only 
sufficient to barely detect a 10 35 erg s _1 magnetar to a 
distance of 4 kpc from Earth. On the other hand , the 
detections of XT E J1810-197 (ICamilo et al.ll2006T l and 
IE 1547.0-5408 (jCamilo et all l2007h in the radio with 
characteristic intensities of >100 mJy kpc 2 suggests that 
future radio surveys with the Low-Frequency Array, the 
Milcu ra Wide-field Array, or the Square K ilometer Array 
(e.g.. ICordes. Lazio. fc McLaughlin! 12004) could provide 
the best constraints on the numbers of magnetars in the 
Galaxy. 

We thank E. Gotthelf for sharing the results of his ob- 
servations of AX J1853. 3-0128. This work made use of 
data obtained from the High Energy Astrophysics Sci- 
ence Archive Research Center (HEASARC), provided by 
NASA's Goddard Space Flight Center. 
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